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PACS 03 . 65 . Wj - State reconstruction, quantum tomography 

Abstract. - We address the intrinsic multimode nature of the quantum state of light obtained 
by pulsed spontaneous parametric downconversion and develop a theoretical model based only 
on experimentally accessible quantities. We exploit the pairwise entanglement as a resource for 
conditional multimode measurements and derive closed formulas for the detection probability 
and the density matrix of the conditional states. We present a set of experiments performed to 
validate our model in different conditions that are in excellent agreement with experimental data. 
Finally, we evaluate nonGaussianity of the conditional states obtained from our source with the 
aim of discussing the effects of the different experimental parameters on the efficacy of this type 
of conditional state preparation. 



Nonclassical states of the radiation field represent a 
resource for quantum information and communication 
and much attention has been devoted to their generation 
schemes. As a matter of fact, beside squeezing, nonclassi- 
cal effects are generally observed in connection with non- 
Gaussian states of light, and this usually implies the pres- 
ence of fluctuating parameters [TJ [2] or of nonlinearities 
higher than the second order are involved (e.g. the Kerr 
effect [3H5]) in the generation scheme. On the other hand, 
the reduction postulate provides an alternative mechanism 
to achieve effective nonlinear dynamics. In fact, if a mea- 
surement is performed on a portion of a composite entan- 
gled system, the other component is conditionally reduced 
according to the outcome of the measurement. The result- 
ing dynamics may be highly nonlinear, and may produce 
quantum states that cannot be generated by currently 
achievable nonlinear processes [6J. The efficiency of the 
process, i.e. the rate of success in getting a certain state, 
is equal to the probability of obtaining a certain outcome 
from the measurement and it may be higher than nonlin- 
ear efficiency, thus making conditional schemes possibly 
convenient even when a corresponding Hamiltonian pro- 



cess exists. 

The nonlinear dynamics induced by conditional mea- 
surements has been analyzed for a large variety of schemes 
[Br TPHrHrET] . including photon addition and subtraction 
schemes [8r H0lfT2l[T3] . optical state truncation of coher- 
ent states [14] , generation of cat-like states [l"5r[P7] , state 
filtering by active cavities [T5][]j|], synthesis of arbitrary 
unitary operators [20 and generation of optical qubit by 
conditional interferometry [21] . More recently, nonGaus- 
sianity of states and operations has been recognized as a 
relevant resource for a series of tasks including entangle- 
ment distillation [28-31] and improvements in both tele- 
portation [8THU 1 I32 ] . cloning [33] and storage [34]. Condi- 
tional state generation has been achieved in the low energy 
regime 35-37] by using single-photon detectors, and the 
question arises whether analogue schemes may be used 
also in the mesoscopic domain [38) . 

In this paper we address multimode conditional mea- 
surements and demonstrate a novel bright source of non- 
classical states [5S] based on (i) pulsed multimode sponta- 
neous parametric down-conversion (PDC) [40TI47] . which 
produces entangled states with a mesoscopic number of 
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photons and (ii) a conditional intensity measurement per- 
formed by photoemissive detectors, called hybrid photode- 
tectors, that are able to partially resolve the number of de- 
tected photons [33 • We develop a theoretical model based 
only on experimentally accessible quantities and derive 
closed formulas for both the detection probabilities and 
the conditional states. We find an excellent agreement 
with the experimental data and succeed in evaluating the 
amount of nonGaussianity of the conditional states despite 
the multimode character of the entangled state. 

The pair of intense correlated beams obtained by pulsed 
PDC represents a convenient system for state prepara- 
tion by conditional measurements. In this case, the state 
outgoing the crystal is intrinsically multimode because 
of the pulsed nature of the pump and the properties of 
the nonlinear interaction [42j;49j, whereas correlations are 
provided by the pairwise entanglement induced by spon- 
taneous PDC. If we assume that the output energy is 
equally distributed among the /i modes of each beam, then 
the overall multimode state produced by pulsed PDC can 
be written as a tensor product of /i identical twin-beam 
states, i.e., 

H = ®|A)) fcfc «A| 
fe=i 



|A)) = V / l-A 2 ^A"|n)®|n) 

n 

with A 2 = N/{n + N), N being the mean total number of 
photons in either of the two beams. In our scheme, which 
is sketched in Fig. [TJ conditional preparation is obtained 
when one of the two beams undergoes a photon counting 
process. If we assume that the detector efficiency 77 is the 
same for each of the /x modes, the probability operator- 
valued measure (POVM){II m } describing the detection 
of m photoelectrons may be written as 

m 

1 3 = 1 

where q = {qi,..., <7 M }, 7 = YX=i *> ^hk is the Kronecker 
delta, and 



k—q ^ ' 

denotes single-mode photon counting POVM. The joint 
probability distribution of photoelectrons is given by 
P\i{s-,t) = Ti'12 [-RII S (g) lit], that, after some algebra, 
reads 
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where M= ±Tr r2 [R E s sII s®Ei fn i] = V N is the total 
mean number of photoelectrons measured on each of the 
two beams. Notice that ([I]) only contains quantities that 
can be experimentally accessed by direct detection. 

When one beam is detected, say the idler, and t pho- 
toelectrons are obtained in the measurement, the cor- 
responding conditional state of the signal is given by 
g t = l/p 2 (t)Tr 2 [RI®n t ] where p 2 (t) = £.Pia(*,*) is 
the marginal probability of measuring t photoelectrons on 
the idler beam. After some calculations we arrive at 
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where 8(x) is the Heaviside step function, 
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and 

M t = T>i[ft E sIIs ] = [tiM+wHvMil-rjWM+fi)- 1 

s 

is the mean number of photoelectrons for the conditional 
state g t . Similarly, upon selecting a set of possible results 
t £ T according to a given rule 7", a suitable engineering of 
the conditional state gj- may be achieved. As an example, 
we will consider the states 

obtained by keeping the photoelectrons on the idler that 
are larger or smaller than a given threshold t*. 
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Fig. 1: (Color online) Schematic diagram of the experimental 
setup. BBOl and BB02: nonlinear crystals; F: cut-off filter; 
HS: harmonic separators; PH: pin-hole apertures; L: lenses; 
MF: multimode optical fibers; HPD: amplified hybrid photode- 
tectors; SGI: synchronous gated integrator; PC: digitizing PC 
board. 

The experimental setup is sketched in Fig. [T] The light 
source was a Nd:YLF ps-pulsed laser (High-Q Laser Pro- 
duction, Austria) with built-in second and third harmonic 
generation. The output at the fundamental (1047 nm) 
and that at the third harmonics (349 nm) were used to 
produce a UV pump field (261.75 nm) via non-collinear 
sum-frequency generation in a BBO crystal (/3-BaB 2 04, 
Castech, China, cut angle 37°, 8 mm length). The pump 
was then sent into another BBO crystal (Kaston, China, 
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cut angle 48.84°, 4 mm length) to produce pairwise en- 
tanglement at 523.5 nm, in order to match the maximum 
quantum efficiency (rj ~ 50%) of our hybrid photodetec- 
tors (Rf 0467U-40, Hamamatsu, Japan). For their photon- 
number resolving power, these detectors proved to be use- 
ful in reconstructing detected-photon distributions [48150] . 
Here we exploit their features to perform both conclusive 
and inconclusive conditional measurements of the photon 
number. The UV stray light was cut-off by a filter and by 
two harmonic separators. Signal and idler were selected 
by two pin-holes (200 or 300 fim diameter, located at 1 m 
from BB02) in order to minimize the number of collected 
modes. Notice that the number of temporal modes, which 
is evaluated from the marginal detected-photon number 
distribution, cannot be reduced at will. The only way to 
reduce the number of modes is to select a single spatial 
mode, which involves the challenging matching of the col- 
lection areas in signal and idler. The possible mismatch 
between the collection areas results in an effective detec- 
tion efficiency, reduced with respect to the nominal effi- 
ciency of the detectors, which can be estimated through 
the level of noise reduction R = <r 2 (s—t)/(s+t) = 1—rj [42] 
exhibited by two beams. For our beams we obtain, with- 
out noise subtraction, 77 ~ 0.06. The light passing the 
pin-holes was coupled to two multimode optical fibers and 
delivered to the detectors, whose outputs were amplified 
(preamplifier A250 plus amplifier A275, Amptek), syn- 
chronously integrated (SGI, SR250, Stanford), digitized 
(ATMIO-16E-I, National Instruments) and, finally, pro- 
cessed off-line. Each experimental run was performed on 
50 000 subsequent laser shots at fixed values of the pump 
intensity. 

As a first test of the correctness of our multimode de- 
scription we checked the expression of pi2(s, t) against 
data: in Fig.[2]we report the experimental joint probabil- 
ity distribution superimposed to the theoretical one, eval- 
uated for the experimental values of the parameters (panel 
(a): PH = 200/jm, /j = 197, 77 = 0.06 and M = 13.4; panel 
(b): PH = 300 /im, fi = 25, r/ = 0.056 and M = 17.1). 
The experimental results fit the theory very well and the 
fidelity ^ st \/p*2 ( s > *)Pi2 P ( s > *) exceeds 0.99 for the whole 
range of parameters. We also notice that the marginal 
probability distributions p\(s) andp2(t) ; are multithermal 
distributions as it has been already observed in experi- 
ments performed at different intensity regimes [42ll48ll49] . 



In Fig.[3]we report the photon distributions Pi\2(s\T) = 
Tri [qt TL g ] of conditional states as obtained from the state 
in Fig. [2] (a) by choosing the values of the measured pho- 
tons on the idler beam according to a given rule and se- 
lecting the corresponding ensemble on the signal beam. 
Panel (a) of Fig. [3] displays the distributions for the de- 
tected photons state Qt obtained by choosing a definite 
number of detected photons (t — 10 and t — 15); panel (b) 
those for gi ' , obtained by keeping the values of detected 
photons larger than a threshold t* (t* = 11 and t* = 17); 





Fig. 2: (Color online) Joint probability distribution of photo- 
electrons pi2(s,t) compared to the experimental points. Panel 
(a): fj, = 197, 77 = 0.06 and M = 13.4; panel (b): fj, = 25, 
77 = 0.056 and M = 17.1). 



finally, panel (c) those for gi~ (t* = 8 and t* = 15). We 
notice that (i) the results are in excellent agreement with 
theory and (ii) despite the small value of effective quan- 
tum efficiency the "conditioning power" of the measure- 
ment (i.e. the differences between the conditional states 
and the corresponding original ones) is appreciable. This 
is clearly illustrated by the behavior of the mean values of 
the distributions, which are reported in panel (d) of Fig. [3] 
as a function of either the conditioning value or the thresh- 
old: the experimental data are in excellent agreement with 
the predictions for M t . 

In Fig. 0] we show the results for another dataset hav- 
ing similar mean value and a considerably lower number 
of modes. We note that the results are again in excel- 
lent agreement with theory and the efficacy of the "condi- 
tioning power" of the measurement is more evident with 
respect to the case of a larger number of modes. 

The nonGaussian character of the conditional states 
may be foreseen from the deviation of the detected-photon 
statistics from that of the original state. However, the 
shape of the distributions (not too different from the un- 
conditioned ones) and the low value of the quantum effi- 
ciency anticipate that the amount of nonGaussianity will 
be unavoidably small. In order to assess the performances 
of our scheme we focus on the conditional state g t and 
use the nonGaussianity measure S[g] — S[t] — S[g], where 
S[g] is the Von Neumann entropy of the state g and r is 
the Gaussian reference state of g, i.e., a Gaussian state 
with the same mean value and covariance matrix as g. S 
has been proved to be a proper measure of nonGaussianity 
[5"T] , as well as a critical parameter to asses nonGaussian- 
ity as a resource [55] . In our case r is a factorized thermal 
state with Mt/ij/j, mean photons per mode [51] and the 
Von Neumann entropy of the conditional state is given 
by S[g t ] = -E"=o ( 7+ ^> t ( 7 )logu; t (7). By using the 
above expression and the Von Neumann entropy of a fac- 
torized thermal state, we evaluate the nonGaussianity and 
normalize its value to that of a maximally nonGaussian 
state for the same mean number of photons and modes, 
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Fig. 3: (Color Online) Photon distribution for conditional 
states, (a): Experimental results (points) and theoretical dis- 
tribution (histograms) for the photoelectrons in the conditional 
signal state gt for t = 10 (red histogram) and t — 15 (gray his- 
togram). The black line and the full circles represent respec- 
tively the theoretical and experimental distribution for the un- 
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conditional state, (b): As in panel (a) for gX^ 1 , t* = 11 
and t* = 17 (gray), (c): As in panel (a) for g* , t* — 8 
(red) and t* — 15 (gray). Panel (d): experimental mean value 
of the distributions as a function of the conditioning value (or 
threshold). Black circles refer to g t , red circles to gl and gray 
circles to g, . Solid lines are the theoretical predictions ob- 
tained for Mt- The dashed line corresponds to the mean value 
of the unconditioned state. The other involved parameters are: 
H = 197, T) = 0.06 and M = 13.4. 
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Fig. 4: (Color Online) Photon distribution for conditional 
states, (a): Experimental results (points) and theoretical dis- 
tribution (histograms) for the photoelectrons in the conditional 
signal state gt for t — 13 (red histogram) and t = 19 (gray his- 
togram). The black line and the full circles represent respec- 
tively the theoretical and experimental distribution for the un- 
conditional state, (b): As in panel (a) for gi , t* — 17 (red) 
and t* = 21 (gray), (c): As in panel (a) for gl~ , t* = 10 
(red) and t* = 15 (gray). Panel (d): experimental mean value 
of the distributions as a function of the conditioning value (or 
threshold). Black circles refer to gt, red circles to gl and gray 
circles to pi . Solid lines are the theoretical predictions ob- 
tained for Mt- The dashed line corresponds to the mean value 
of the unconditioned state. The other involved parameters are: 
H = 25, n = 0.056 and M = 17.1. 



i.e., a factorized Fock state [5T]. The renormalized non- 
Gaussianity Sn[g t ) is reported in Fig. [5]for different values 
of the experimental parameters. 

As it is apparent from the plots the renormalized non- 
Gaussianity Sr is a decreasing function of the energy of the 
conditional state and of the number of modes. The effect 
of the quantum efficiency is more relevant for large number 
of modes and large energy. As it concerns the conditioning 
value t of detected photons, we see that Sr monotonically 
increases with t, again with the quantum efficiency playing 
a major role for large number of modes. Overall, the goal 
of achieving high nonGaussianity requires a small num- 
ber of modes or, at fixed number of modes, a high value 
of the quantum efficiency. Since the pulsed nature of the 
PDC pump unavoidably leads to a multimode output, the 
performances of the present source in the generation of 
nonGaussian states may be improved by increasing the 
overall quantum efficiency, i.e. the matching of the collec- 
tion areas in signal and idler beams. 

In conclusion, we have suggested and demonstrated a 
novel bright source of nonclassical states based on mul- 
timode spontaneous PDC and conditional intensity mea- 
surements. We have developed a theoretical model based 
only on experimentally accessible quantities and derived 
closed formulas for the detection probability, the condi- 
tional states and the corresponding nonGaussianity. We 



have compared our predictions with experimental data 
and found an excellent agreement in the whole range of 
accessible experimental parameters. Our results clearly 
indicate the possibility of quantum state engineering with 
multiphoton/multimode conditional states using meso- 
scopic photon counting and multimode pairwise correlated 
states. 

This work has been partially supported by the CNR- 
CNISM agreement. 



REFERENCES 

[1] V. D'Auria, C. de Lisio, A. Porzio, S. Solimeno, J. Anwar 

and M. G. A. Paris, Phys. Rev. A, 81 (2010) 033846. 
[2] A. Chiummo, M. De Laurentis, A. Porzio, S. Solimeno 

and M. G. A. Paris, Opt. Expr., 13 (2005) 948. 
[3] C. Silberhorn, P. K. Lam, O. Weifi, F. Konig, N. Ko- 

rolkova, and G. Leuchs, Phys. Rev. Lett, 86 (2001) 4267. 
[4] O. Glockl, U. L. Andersen and G. Leuchs, Phys. Rev. A 

73 (2006) 012306. 
[5] T. Tyc and N. Korolkova, New J. Phys., 10, (2008) 

023041. 
[6] M. G. A. Paris, M. Cola and R. Bonifacio, Phys. Rev. A, 

67, (2003) 042104. 
[7] G. M. DAriano, P. Kumar, C. Macchiavello, L. Maccone 

and N. Sterpi, Phys. Rev. Lett, 83 (1999) 2490. 



p-4 



Conditional measurements on multimode pairwise entangled states from spontaneous parametric downconversion 




Fig. 5: (Color Online) Renormalized nonGaussianity Sr of the 
conditional state Qt for different values of the experimental pa- 
rameters. (Top left panel): Sr as a function of the total num- [25 
ber of measured photoelectrons M t for t = 5 and for different [26 
values of the quantum efficiency rj and of the number of modes 
(i. Solid lines are for ji = f 97, dotted for it, — 25 and dashed for [27 
ji, — 1 . For each group of lines we have, from bottom to top, [28 
rj = 6% (red), r) = 8% (green), rj = 10% (blue) and rj = 20% 
(black). (Top right panel): Sr as a function of the condition- [29 
ing number t of detected photoelectrons on the idler beam for [30 
Mt = 4 and for different values of the quantum efficiency rj and 
of the number of modes ji. Solid lines are for jjl = 197, dot- 
ted for fj, = 25 and dashed for ii — 1. For each group of lines [31 
we have, from bottom to top, rj — 6% (red), rj = 8% (green), 
rj = 10% (blue) and r) = 20% (black). ( Bottom left panel ): Sr 
as a function of the quantum efficiency rj for M t — 4 and for [32 
different values of the conditioning number t of detected pho- 
toelectrons and of the number of modes fj,. Solid lines are for [33 
t — 5 and dashed ones for i = 15. From bottom to top we can 
see lines for fi — 197 (blue), fi = 25 (green), and ji — 1 (red). [34 
(Bottom right panel): Sr as a function of the number of modes 
fi for Mt = 4 and for different values of conditioning number [35 
t of detected photoelectrons and for the quantum efficiency rj. 
Solid lines are for t = 2 and dashed for t = 15. For each group [36 
of lines we have, from bottom to top, rj = 6% (red), rj = 8% 
(green), rj = 10% (blue) and rj = 20% (black). 

[37 

[8] T. Opatrny, G. Kurizki and D. -G. Welsch, Phys. Rev. 

A, 61, (2000) 032302; M. Dakna, L. Knoll and D. -G. [38; 

Welsch, Opt. Comm., 145 (1998) 309. [39 

[9] P. T. Cochrane, T. C. Ralph and G. J. Milburn, Phys. 

Rev. A, 65 (2002) 062306. 
[10] S. Olivares, M. G. A. Paris and R. Bonifacio, Phys. Rev. [40 

A, 67 (2003) 032314; S. Olivares, M. G. A. Paris, J. Opt. 

B, 7 (2005) S392. [41 
[11] S. Olivares and M. G. A. Paris, J. Opt. B, 7 (2005) 616; 

C, Invernizzi, S. Olivares, M. G. A. Paris and K. Ba- [42 
naszek, Phys. Rev. A, 72 (2005) 042105. 

[12] P. Marek and R. Filip, Phys. Rev. A, 81 (2010) 022302. [43; 
[13] A. Allevi, A. Andreoni, M. Bondani, M. G. Genoni, S. 



Olivares, Phys. Rev. A (2010), in press. 

D. T. Pegg, L. S. Philips and S. M. Barnett, Phys. Rev. 

Lett, 81 (1998) 1604. 

M. Dakna, T. Anhut, T. Opatrny, L. Knoll and 

D. G. Welsch, Phys. Rev. A, 55 (1997) 3184. 

M. Dakna, J. Clausen, L. Knoll and D.-G. Welsch, Acta 

Phys. Slov., 48 (1998) 207. 

S. B. Zheng, Phys. Lett. A, 245 (1998) 11. 

G. M. D'Ariano, L. Maccone, M. G. A. Paris and 

M. F. Sacchi, Phys. Rev. A, 61 (2000) 053817; Fort. Phys., 

48 (2000) 511. 

Lu-Ming Duan, G. Giedke, J. I. Cirac and P. Zoller, Phys. 

Rev. Lett, 84 (2000) 4002. 

B. Hladky, G. Drobny and V. Buzek, Phys. Rev. A, 61, 

(2000) 022102. 

M. G. A. Paris, Phys. Rev. A, 62 (2000) 033813. 

J. Clausen, M. Dakna, L. Knoll and D.-G. Welsch J. Opt. 

B, 1 (1999) 332; M. G. A. Paris, Phys. Lett. A 217 (1996) 
78. 

A. Napoli, A. Messina and S. Maniscalco, Acta Phys. 
Slov., 50, (2000) 519. 

F. Plastina and F. Piperno, Eur. Phys. J. D, 5 (1999) 
411. 

M. Ban, Opt Comm., 143 (1997) 225. 

A. Kozhekin, G. Kurizky and B. Sherman, Phys. Rev. A, 

54, (1996) 3535. 

M. S. Kim, J. Phys. B, 41 (2008) 133001. 

J. Eisert, S. Scheel and M. B. Plenio, Phys. Rev. Lett. 89 

(2002) 137903. 

J. Fiurasek, Phys. Rev. Lett, 89 (2002) 137904. 

R. Dong, M. Lassen, J. Heersink, C. Marquardt, R. Filip, 

G. Leuchs and U. L. Andersen, Nature Phys., 4 (2008) 
919. 

T. Aoki, G. Takahashi, T. Kajiya, J. Yoshikawa, 
S. L. Braunstein, P. van Loock and A. Furusawa, Nature 
Phys., 5 (2009) 541. 

C. Invernizzi, S. Olivares, M. G. A. Paris and K. Ba- 
naszek, Phys. Rev. A, 72 (2005) 042105. 

N. J. Cerf, O. Kriiger, P. Navez, R. F. Werner and 

M. M. Wolf, Phys. Rev. Lett, 95 (2005) 070501. 

F. Casagrande, A. Lulli and M. G. A. Paris, Phys. Rev. 

A, 75, (2007) 032336. 

A. Ourjoumtsev, A. Dantan, R. Tualle-Brouri and 

P. Grangier, Phys. Rev. Lett, 98 (2007) 030502. 

H. Takahashi, J. S. Neergaard-Nielsen, M. Takeuchi, 

M. Takeoka, K. Hayasaka, A. Furusawa and M. Sasaki, 

Nature Phot, 4 (2010) 178. 

T. Gerrits, S. Glancy, T. S. Clement, B. Calkins, 

A. E. Lita, A. J. Miller, A. L. Migdall, S. W. Nam, 

R. P. Mirin and E. Knill, arXive:1004.2727vl 

M. Sasaki and S. Suzuki, Phys. Rev. A, 73 (2006) 043807. 

M. G. Genoni, F. A. Beduini, A. Allevi, M. Bondani, 

S. Olivares and M. G. A. Paris, Phys. Scripta, T139 

(2010) in press 

A. Agliati, M. Bondani, A. Andreoni, G. De Cillis and 

M. G. A. Paris, J. Opt. B, 7,( 2005) 652. 

W. Wasilewski, A. I. Lvovsky, K. Banaszek, C. Radzewicz. 

Phys. Rev. A 73, (2006) 063819. 

M. Bondani, A. Allevi, G. Zambra, M. G. A. Paris and 

A. Andreoni, Phys. Rev. A, 76 (2007) 013833 . 

A. I. Lvovsky, W. Wasilewski, K Banaszek, J. Mod. Opt 

54, (2007), 721. 



p-5 



A. Allevi et al. 



[44; 
[45; 

[4e; 
[47; 

[48 
[49 
[50 
[51 

[52 



W. Wasilewski, C. Radzewicz, R. Frankowski, K. Ba- 
naszek, Phys. Rev. A, 78 (2008) 033831. 
W. Mauerer, M. Avenhaus, W. Helwig and C. Silberhorn, 
Phys. Rev. A, 80 (2009) 053815. 

W. Helwig, W. Mauerer and C. Silberhorn Phys. Rev. A 
80, (2009) 052326. 

C. Soller, B. Brecht, R J. Mosley, L. Y. Zang, 
A. Podlipensky, N. Y. Joly, P. St. J. Russell and C. Sil- 
berhorn, Phys. Rev. A, 81, (2010), 031801. 
M. Bondani, A. Allevi, A. Agliati and A. Andreoni, J. 
Mod. Opt. 56 (2009) 226. 

F. Paleari, A. Andreoni, G. Zambra and M. Bondani, Opt. 
Express 12 (2004) 2816. 

M. Bondani, A. Allevi and A. Andreoni, Adv. Set. Lett., 
2 (2009) 463; Opt. Letters, 34 (2009) 1444. 
M. G. Genoni, M. G. A. Paris and K. Banaszek, Phys. 
Rev. A 76 (2007) 042327; Phys. Rev. A, 78 (2008) 
060303(R). 

M. G. Genoni, C. Invernizzi and M. G. A. Paris, Phys. 
Rev. A, 80 (2009) 033842. 



p-6 




x10 
10 

8 

E 6 



-2 



(c) 



<s. 



4 



§ 




if 


* 


J P 


f 


tf 





» 

D 









o I 



d 
1 




10 15 
s 











rill I II h>;; 



5 



10 



x10 



-2 



15 



r 10 



T* 

Q. 



5 











\0 

I 



•o 



5 



t 
t 




